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ABSTRACT 

We present results of near-IR spectroscopic measurements of 7 star-forming galaxies at 2.1 < z < 2.5. 
Drawn from a large spectroscopic survey of galaxies photometrically pre-selected by their U n GlZ colors 
to lie at z ~ 2, these galaxies were chosen for their bright rest-frame optical luminosities (K s < 20.0). 
Most strikingly, the majority of the sample of 7 galaxies exhibit [Nil] / Ha nebular emission line ratios 
indicative of at least solar H II region metallicities, at a lookback time of 10.5 Gyr. The broadband 
colors of the iV-bright sample indicate that most have been forming stars for more than a Gyr at z ~ 2, 
and have already formed stellar masses in excess of 10 11 M Q . The descendants of these galaxies in 
the local universe are most likely metal-rich and massive spiral and elliptical galaxies, while plausible 
progenitors for them can be found among the population of z ~ 3 Lyman Break Galaxies. While 
the iis-bright z ~ 2 galaxies appear to be highly evolved systems, their large Ha luminosities and 
uncorrected Ha star- formation rates of 24 — 60A/ Q yr~ 1 indicate that active star formation is still 
ongoing. The luminous UV-selected objects presented here comprise more than half of the high- 
redshift [z > 1.5) tails of current iVband-selected samples such as the K20 and Gemini Deep Deep 
surveys. 

Subject headings: galaxies: evolution — galaxies: high-redshift - galaxies: abundances 



1. INTRODUCTION 

Rest-frame optical emission lines from extragalactic 
H II regions provide an important diagnostic of galac- 
tic chemical evolution and therefore serve as an inde- 
pendent gauge of the history of star formation in the 
universe. Furthermore, the chemical evolution of galax- 
ies and of the IGM are interconnected processes since 
feedback from star formation in galaxies may be re- 
sponsible for polluting the IGM with heavy elements 
over a wide range of redshifts ijHeckman et alJ l200Ct 
lAdelberger et all 12003). Now, with the advent of sen- 
sitive near-infrared spectrographs on 8 — 10 m class tele- 
scopes, it is possible to measure rest-frame optical neb- 
ular emission lines out to z > 3, and therefore to trace 
galactic chemical evolution over the bulk of cosmic his- 
tory. 

In the local universe, star-forming spiral and irregu- 
lar galaxies display a strong correlation between their 

1 Based, in part, on data obtained at the W.M. Keck Observa- 
tory, which is operated as a scientific partnership among the Cal- 
ifornia Institute of Technology, the University of California, and 
NASA, and was made possible by the generous financial support 
of the W.M. Keck Foundation. 



H II region oxygen abundances and optical luminosi- 
ties, according to which more luminous and massive 
galaxies exhibit a highe r degree of metal enrichment 
(iGarnett fc Shields! 119871 iGarnettJ 120021: iZaritskv et alJ 
119941 lLamareille et al.ll2004t iTremonti et al.ll2004ft . An 
analogous correlation is observed in elliptical galax- 
ies, where chemical abundances are measured from the 
strengt hs of stellar magnesium and iron absorption line 
indices l)Brodie fc Huchr g||l991). As we move to higher 
redshifts, the metallicity-luminosity relationship ma y 
evolve in both zero point and slope (jMaier et alj f20041. 
Using a sample of 64 star-forming galaxies at 0.26 < 
z < 0.82 drawn from the Deep Extrag alactic Evolution- 
ary Pro be Groth Strip Survey fDGSSI. lKobulnickv et all 
(2003) demonstrate that galaxies at 0.6 < z < 0.8 
are 1 — 3 magnitudes more luminous than local galax- 
ies of similar metallicity. In contrast. iLillv et alJ (2003) 
find that a majority of a sample of 66 star-forming 
galaxies at 0.47 < z < 0.92 drawn from the Canada- 
France Redshift Survey (CFRS) have similar metallic- 
ities to galaxies in the local universe with comparable 
luminosities. At much higher redshift, z > 2.5, most 
of the information about the abundances of galaxies 
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has been derived from absorption li nes in Damped Ly- 
man Alpha systems (DL As) (e.g., iPettini 
iProchaska Sz WoTfel l2002j) . However, relating line-of- 
sight DLA absorption abundances to abundances mea- 
sured from global galactic H 11 emission-line spectra is 
not trivial, at least in part because the relationship of 
DLAs themselves to high-redshift UV-selected galaxies 
has not been established. There are only limited data 
on the H 11 region oxygen abundances of star-forming 
galaxies at z > 2.5. A sample of six galaxies has been as- 
sembled with near-infrared spectroscopic observations of 
the strong nebular emission lines that enable abundanc e 
estimates llPettini et a,l .1127)011 IKohulnickv fc Kooll200fft. 
While there are only rough constraints on the abun- 
dances in these galaxies, with (O/H) ~ 0.1 — 1.0 (O/H) , 
they are seen to be 2 — 4 magnitudes more luminous than 
local galaxies of similar metallicity. 

Until now, there have been virtually no observations of 
the H II region metallicities in star-forming galaxies at 
1.4 < z < 2.5, even though this critical epoch may host 
the production of a large fraction of the heavy elements 
present in the local universe. Galaxies in the "redshift 
desert" range have only recently become accessible for 
detailed study due, in large part, to improvements in 
both near-UV and 0.8 — 1.0/im spectroscopic capabilities 
on 8—10 m class telescopes. Employing new technologies, 
several current surveys utilizing complementary selection 
techniques contain spectroscopically confirmed galaxies 
at thes e redshifts. The sel e ction criteria include op tical 
colors (IRteidel et all l200l lAdelberger et all l200l. K- 
band magnitude ( Cim atti et al.ll2002t iGlazcbroo k et all 
|2004(), J—K color llFranx et al.l2003j) . and submillimetcr 
flux IChapma n et al.ll20 03K A complete picture of the 
galaxy population at z ~ 2 will require an understand- 
ing of how these several selection techniques overlap and 
complement ea ch other. 

As s hown bv IPettini fc Pagell (|2004Tl and lSteidel et all 
( 2004) , H II region metallicities can be measured in red- 
shift desert galaxies using the ratio of the [Nil] A6584 
and Ha nebular emission lines. Here we present near- 
infrared (i-T-band) spectroscopic observations of [Nil] 
and Ha emission for a small sample of UV-selected star- 
forming galaxies at 2.1 < z < 2.5. These galaxies were 
selected for observation because they are luminous in 
the rest-frame optical (K s < 20.0), and therefore ap- 
plication of the same observational strategy yields much 
higher signal-to-noise (S/N) spectra than previous z ~ 2 
nea r-infrared spectr oscopic observations of fainter galax- 
ies (|Erb et al.1 12003(1 . Along with the spectroscopic ob- 
servations, which probe the metal content of the galaxies, 
we also use broad-band optical/near- infrared colors to es- 
timate their stellar populations and masses. With such 
information, it is possible to consider the build-up of stel- 
lar mass in the early universe, and the epoch at which 
the progenitors of massive galaxies in the local universe 
were mostly assembled. 

The sample of z ~ 2 galaxies and near-infrared spec- 
troscopic observations are described in |J21 In fJ3 the re- 
sults of the Ha and [Nil] spectroscopic measurements 
are presented. Oxygen abundances are derived from 
these results in jJU In the discussion is extended to 
the galaxy metallicity-luminosity relationship, the stel- 
lar populations and masses of the luminous z ~ 2 sam- 
ple, a comparison of this sample with other recent sur- 
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Fig. 1. — The Color-magnitude diagram of 121 UV-selected 
galaxies at z > 1 ((«) = 2.28 ± 0.28), in the Q1623 field. Tl 
and K s magnitudes are on the AB and Vega systems, respectively. 
The dotted diagonal line indicates the survey magnitude limit of 
■R = 25.5. Only ~ 10% of the objects at z > 1 also have K s < 20.0, 
indicated by the long-dashed line. The nine K < 20.0 objects 
observed with NIRSPEC in September 2003 are indicated by large 
blue triangles. The object in the lower left-hand corner, with K a = 
18.29 and U- K a = 2.15, is a QSO at z = 2.529. The median color 
uncertainty is cr(TZ — K s ) = 0.25, though the 10 _R" 3 -selected objects 
have a{Tl - K s ) ~ 0.10. 



veys, and a comment on inferring stellar masses from 
population synthesis modeling versus dynamical masses 
from Ha line- widths. Finally, in fj^J we summarize our 
main conclusions. For all luminosity, star-formation rate, 
and stellar mass calculations, we adopt a cosmology with 
H = 70 km s- 1 Mpc -1 , n m = 0.3, and A = 0.7. 

2. SAMPLE AND OBSERVATIONS 

The objects di scussed here were d rawn from the sur- 
vey presented bv lSteidel et alJ (J.2004). This survey pho- 
tometrically identifies star-forming galaxies at z ~ 2 
by their locations in U n — G vs. G — 1Z color space. 
Objects satisfying the o ptical color criteria presented in 
lAdelberger et al] l)2004|) have been followed up spectro- 
scopically using the Low Res olution Imaging Spectrom- 
eter on the Keck I telescope (|Oke et al.l ll995). To date, 
692 galaxies have been spectroscopically confirmed with 
1.4 < z < 2.5. The survey has been conducted in 
seven separate fields, five of which contain one or more 
bright background QSOs at z }t 2.5 in order to study 
the evolving relationship between star- forming galaxies, 
gas and metals in the intergalactic medium (IGM). The 
remaining two pointings are the GOODS-N and West- 
phal/Groth fields, both of which have extensive multi- 
wavelength datasets that are or will soon be publicly 
available. In addition to optical imaging and spectro- 
scopic data, in June and October 2003 we obtained near- 
infrared i^T s -band photometry for galaxies in three of the 
survey QSO fields using the Wide Field Infrared Camera 
(WIRC) on the Palomar 5.1m Hale telescope. K s magni- 
tudes have been measured for 283 spectroscopically con- 
firmed galaxies at z > 1.4, with (z) = 2.25±0.31. One of 
these three fields, Q1623 (1623+27), contains 167 spec- 
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Table 1. Galaxies Observed With Keck II/NIRSPEC 



Galaxy 


R. 


A. 


(J2000) 


Dec. 


(J2000) 






ZHc C 


TIab 


(G - TV)ab 


(U n — G)ab 


TlAB — Ks.Vega 


Exposure (s) 


Q1623-BX274 


16 


25 


38.202 


26 45 


57.141 


2.415 


2.408 


2.4100 


23.23 


0.25 


0.89 


3.66 


3x900 


Q1623-MD66 


1(5 


25 


40.392 


26 50 


08.878 




2.111 


2.1075 


23.95 


0.37 


1.40 


4.04 


3x900 


Q1623-BX341 


16 


25 


43.554 


26 46 


36.942 




2.377 




24.83 


0.47 


0.90 


4.81 


2x900 


Q1623-BX344 d 


16 


25 


43.931 


26 43 


41.977 




2.422 


2.4224 


24.42 


0.39 


1.25 


4.41 


2x900 


Q1623-BX453 


16 


25 


50.836 


26 49 


31.399 


2.183 


2.171 


2.1816 


23.38 


0.48 


0.99 


3.65 


3x900 


Q1623-BX513 


1(5 


25 


55.856 


26 46 


50.304 


2.249 


2.244 


2.2473 


23.25 


0.26 


0.68 


3.33 


2x900 


Q1623-BX528 


1(5 


25 


56.439 


26 50 


15.444 




2.266 


2.2682 


23.56 


0.25 


0.71 


3.87 


4x900 


Q1623-BX599 


16 


2(5 


02.545 


26 45 


31.900 




2.329 


2.3304 


23.44 


0.22 


0.80 


3.65 


4x900 


Q1623-BX663 C 


1(5 


26 


04.576 


26 48 


00.202 


2.435 




2.4333 


24.14 


0.24 


1.02 


4.26 


3x900 



a Vacuum heliocentric redshift of Lya emission line, when present. 

b Vacuum heliocentric redshift from rest-frame UV interstellar absorption lines. 

c Vacuum heliocentric redshift of Ho emission line. 

d BX344 was observed with a 0'.'57 slit, while all other galaxies were observed with a 0'.'76 slit. 

e The Ho emission in BX663 has a two distinct peaks, one at z = 2.4333, as listed above, and another at z = 2.4289. 



troscopically confirmed galaxies at z > 1.4 in the region 
covered by WIRC, 121 of which were detected at K s . 

At the redshifts probed in the sample, if s -band cor- 
responds most closely to rest-frame optical 7?.-band, and 
for z > 1.9 i^-band spectroscopic observations cover the 
Ha and [N II] nebul ar emission lines. Using NIRSPEC 
ijMcLean et al.llT998|) on the Keck II telescope, we have 
assembled a sa mple of more than 50 Ha measurements 
in the K s -b&nd <|Erb et alJl2003ll2004j) . Since the major- 
ity of these objects were chosen pr imarily for their close 
proximity to a QSO sightline (see lErb et al.l l2003. for a 
discussion), the NIRSPEC sample contains a fairly un- 
biased view of the range of rest-frame UV and optical 
properties of galaxies in the lSteidel et alJ l|2004j) survey. 
With this unbiased sample (with respect to photomet- 
ric properties) in place as a control, we used the K s 
photometry to select galaxies for NIRSPEC rest-frame 
optical spectroscopic follow-up. In Figure ^ we present 
the distribution of K s magnitudes and 7Z — K s colors 
for spectroscopically confirmed galaxies in the Q1623 
field. As shown in the figure, 10% of the galaxies in 
the Q1623 field with WIRC photometry h ave K s < 20.0 
(and would be included in the K20 survey ijCimatti et alJ 
2002), the high-rcdshift tail of which overlaps with the 
UV-selected sample here under consideration) , while 8% 
have H — K 8 > 4.0. In order to measure the rest-frame 
optical spectra of objects with bright rest-frame optical 
luminosities, we targeted galaxies with K s < 20.0. 

The Keck II/NIRSPEC observations were conducted 
during the course of a single observing run spanning 9 
- 13 September 2003. We targeted a total of nine spec- 
troscopically confirmed z ~ 2 objects in the Q1623 field 
satisfying K s < 20.0, four of which also had 1Z—K S > 4.0. 
All but one of th ese galaxies satisfy the " BX" color cri- 
teria presented in lAdelberger et ail l|2004fl . while the re- 
maining galaxy falls in the adjacent region of U n GlZ colo r 
space occupied by "MD" objects JStcidel et al. 2003). 
The single MD object in the sample falls on the border 
between the two color criteria, lies at the same redshift as 
the other eight galaxies, and is therefore treated identi- 
cally for the remainder of the analysis and interpretation. 
We followe d the same observ ational procedures described 
in detail in lErb et all ((2003). In brief, we used the NIR- 
SPEC6 filter, which transmits from 1.84 — 2.63 /im. Since 



the medium-dispersion mode in which we were observ- 
ing only allows for ~ 0.4 /im of coverage to fall on the 
NIRSPEC detector at once, we set the NIRSPEC cross- 
disperser angle to cover 1.9 — 2.3 /im, which encompassed 
Ha and [Nil] for the redshifts of all of our target galaxies. 
For all objects except BX344, we used a 0'.'76 x 42" long- 
slit, affording a spectral resolution of ~ 15 A, whereas 
for BX344, we used a 0"56 x 42" long-slit, affording 
~ 12 A spectral resolution, as measured from the widths 
of skylines. 2 The same method of blind-offset objec t ac- 
quisition was used as described in lErb et alJ i|2003f) . In 
most cases, the slit position angle was determined by 
the attempt to fit two galaxies on the slit (usually our 
primary near-IR-selected target and a nearby compan- 
ion with or without a Keck I/LRIS-confirmed redshift), 
though in the case of BX513, which appeared elongated 
in the WIRC K s image, the slit position angle was chosen 
to lie along the major axis of the galaxy light. Targets 
were observed for 2, 3, or 4 x 900 seconds. We detected 
eight out of nine of our target galaxies, and the one unde- 
tected galaxy, BX341, was only observed for 2 x 900 sec- 
onds. Conditions were photometric throughout the run 
with approximately 0'.'5 FWHM seeing in the NIRSPEC6 
band. A summary of the observations including target 
coordinates, redshifts, optical and near-IR photometry, 
and total exposure times is given in Table ^ The two- 
dimensional galaxy spectral images were then reduced, 
extracted to one dimension along with X-a error spec- 
tra, and flux-calibrated with observations of A stars, ac- 
cordin g to the procedures outlined in detail in lErb et alJ 
lf20f)l . 

3. RESULTS 

The one-dimensional, flux-calibrated NIRSPEC spec- 
tra are shown in Figure El Ha and [Nil] A 6584 emis- 
sion are very significantly detected in all spectra except 

2 BX344 was the first near-IR-selected object targeted during the 
observing run, and we used it to experiment with a narrower slit 
setup. A narrower slit reduces the effects of bright sky lines, both 
due to the reduced sky flux in the lines themselves, and the reduced 
number of wavelengths affected by sky flux because of increased 
spectral resolution. However, we found that difficulties in centering 
the object on the narrower slit and the resulting variable flux losses 
outweighed the advantages of the narrower slit. Therefore, the 
remainder of the near-IR selected objects were observed with the 
0'.'76 slit. 
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that of BX513. The redshift of this object places the Ha 
emission line directly on top of a sky OH emission feature 
at 2.1318 /zm, thus making flux measurements extremely 
uncertain. We therefore exclude BX513 from all subse- 
quent analysis, leaving us with a robust sample of seven 
galaxies. 3 [SII] emission is also significantly detected in 
the spectra of MD66 and BX453, and marginally de- 
tected for BX274, BX344, BX599, and BX528. Ha and 
[Nil] emission line fluxes were determined by first fitting 
a one dimensional Gaussian profile to the higher S/N 
Ha feature to obtain the redshift and FWHM (in wave- 
length) . The Ha redshift and FWHM were then used to 
constrain the fit to the [Nil] emission line. This proce- 
dure is based on the assumption that the Ha and [Nil] 
lines have exactly the same redshift and FWHM, but 
that the Ha line offers a higher S/N estimate of these 
parameters. Such an assumption appears to describe all 
of the spectra well except that of BX453, in which the 
[Nil] line has broader wings than the Ha line, and the 
fit based on the Ha parameters u nderestimates the [Nil] 
flux by ~ 30%. lErb et all l|2004D find an average veloc- 
ity dispersion of (a) = 114 km s _1 for the sample of 
61 UV-selected z ~ 2 galaxies with Ha measurements 
from which the current sample is drawn. Five of the 
seven if s -selected objects have a > 114 km s _1 , though 
two of them, BX453 and B X344, have veloci ty disper- 
sions significantly below the lErbeT al. (2004) average. 
We list the measured Ha velocity dispersions, Ha and 
[Nil] fluxes, [Nil]/ Ha flux ratios, and random statistical 
uncertainties on these quantities in Table |3 

Even with relatively short exposure times, we ob- 
tained highly significant detections of Ha emission for 
our target objects. The mean Ha flux is 1.34 x 
10~ 16 erg s _1 cm -2 with a standard deviation among 
the seven galaxies of 0.50 x 10~ 16 erg s _1 cm -2 ; the 
corresponding average Ha star-for mation rate i s 47 ± 
15MQyr -1 using the calibration of iKennicuttl |l998). 
These values are almost three times larger than the av- 
erage Ha flu x and sta r -form ation rate of the sample 
considered in lErb et all l)2003[) . The average rest-frame 
optical continuum luminosity for these galax ies is also 
rough ly three times brighter than that of the lErb et aTl 
(2003) sample, so the average Ha equivalents widths of 
the two samples are very similar. Even after correcting 
the broadband K s magnitudes for the contribution from 
Ha emission lines, five of seven galaxies in our sample 
would still have been selected as K 8 ^ 20, whereas the 
remaining two galaxies, BX344 and MD66, would have 
almost made the cut, with K s = 20.20 and K s — 20.08, 
respectively. In other words, for five of the seven galaxies 
in the sample presented here, the Ha emission line con- 
tributes less than 0.15 of the flux in the if s -band. For the 
remaining two galaxies, the Ha line contributes between 
0.15 and 0.20 of the if s -band flux. 

Most of the galaxies have spatially unresolved Ha pro- 
files, with the exception of BX528 and BX663. As shown 
in Figure|31 BX528 has a spatially resolved and tilted Ha 
line, even though the slit position angle was not selected 
to lie along any morphologically- preferred axis. The ex- 
tent of the velocity shear in BX528 is v max — v m i n ~ 

3 While [Nil] is not detected in the BX513 spectrum, the limits 
on [Nil]/ Ho are not very restrictive, given the low S/N and sky 
contamination of the Ha line. 



140 km s _1 . BX663 has a double-peaked Ha emission 
line, with a velocity separation of Av — 390 km s _1 be- 
tween the two peaks. The higher-redshift Ha component 
is ~ 3 times more luminous than the lower-redshift com- 
ponent and fairly extended, with detected flux extending 
~ 2" along the slit. While there is a tentative detection 
of [Nil] for the fainter, lower-redshift Ha component, 
we only list the [Nil] / Ha ratio for the more-significantly 
detected, higher-redshift component. The Ha flux for 
BX663 listed in Table [21 represents the sum of the two 
components. 

The most striking result, shown in Figure and Ta- 
ble 13 is the high [Nil]/ Ha ratios measured for the K s - 
selected galaxies. The average ratio is [Nil] /Ha = 0.27, 
and the measured values range from 0.17 — 0.43. It 
is relevant to consider the lowest [Nil] / Ha ratio that 
could hav e been measured , given an Ha S/N. In the 
sample of lErb et all l)2003j) . while [Nil] is detected for 
only one object (Q1700-BX691, which has an above- 
average Ha flux), the [Nil] non-detections for the re- 
mainder of the sample are not very restrictive on the 
value of [Nil] /H a. Since the Ha emission lines in the 
lErb et all l)2003t) sample are typically detected with a sig- 
nificance of ~ 10 times the random noise, [Nil] will only 
be significantly detected if it attains ~ 0.30 the strength 
of Ha. To estimate the typical [Nil] / Ha for the galaxies 
in Er b et all l)2003l) . we constructed an average compos- 
ite spectrum from the individual NIRSPEC spectra, and 
found ([Nil]/ Ha) = 0.10. In the ^-selected sample, 
the higher significance of the Ha detections allows for 
a more sensitive probe of individual galaxies' [Nil]/ Ha 
ratios. We could have detected [Nil]/ Ha < 0.10 for 
BX453, MD66, and BX599, and [Nil]/ Ha = 0.12 - 0.13 
for BX274, BX344, and BX663. BX528 is the only galaxy 
whose [Nil]/ Ha ratio of 0.19 is at the limit of what we 
could have detected. 

4. THE OXYGEN ABUNDANCE 

The most robust measurements of chemical abun- 
dances from nebular emission lines are derived 
from the temperature-sensitive ratio of an upper- 
and intermediate-level transition from the same ion. 
The ratio of transitions such as [OIII] A4363 to 
[OIII] AA5007, 4959, indicates the nebular electron tem- 
perature (T e ), and, along with the electron density, yields 
a direct estimate of the chemical abundance. Unfortu- 
nately, the relatively weak [OIII] A4363 line is only ob- 
served in the hottest, most metal-poor H II regions, and 
is not detected in more metal-rich environments. To ad- 
dress this issue, H II region photoionization models that 
predict the relative strengths of Ha, H/3, and strong 
forbidden lines of [Oil], [OIII], [Nil], [SII], and [Sin], 
have been calibrated against H II regions with direct 
abundance estimates llPagel et al.lll979UEvans fc Dopital 
119851 IKewlev fc Dopital l2002|) "~With these empirically 
calibrated models, it is possible to infer approximate 
chemical abundances over a wide range of metallicity, 
using only the ratios of strong, easily detectable H II 
region emission lines. 

The most commonly-used strong-line abundance indi- 
cator is the "-R23" method, which relates the nebular 
oxygen abundance to the ratio of ([Oil] + [OIII])/ H/3 
()Pagel et al.l I1979J) . The R 23 method has been em- 
ployed to determine the metallicities of individual H II 
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Fig. 2. — Fully-reduced one-dimensional NIRSPEC spectra. The Ha, [Nil], and [SII] lines are marked in each spectrum. The solid line 
indicates flux, while the dotted line represents the error spectrum. The vertical scale indicates flux-density in units of 10 — 18 erg s — X A . 
We exclude BX513 (upper right hand panel) from all analysis, since its Ho emission line falls directly on top of a sky line. 



regions in the local universe, and the average metallici- 
ties of entire star-forming galaxies in the distant universe 
Ipettini et alJ2001l:lLiriv et al.l200^ . The double- valued 
nature of the R23 calibration limits its utility unless an 
independent line-ratio can be used to determine whether 
a galaxy is on the upper, metal-rich branch or the lower, 
metal-poor branch. Furthermore, in objects for which 
[OIII] A4363 or other temperature-sensitive auroral lines 
arc detected, mctallicities inferred from R23 are system- 
atically higher by 0.2 — 0.5 dex than those comp uted 
directly from the T e method ijKennicutt et alJfeOOSjh 

As presented in section we have measured [Nil] / Ha 
from the sample of UV-selected star- forming z ~ 2 galax- 
ies. This line ratio can also be used to infer approxi- 
mate chemical abundances. Unlike R23, the [Nil]/ Ha 



ratio increases monotonically with increasing metallic- 
ity (at least up to solar metallicity) . Another benefit 
of the [Nil]/ Ha ratio is the small difference in wave- 
length between the two features, which reduces the im- 
portance of accurate flux-calibration and dust extinc- 
tion correction. The utility of the "iV2" index (N2 = 
log([NII]A 6584/ Ha)) has r ecentl y been revisited empir- 
ically by iPettini fc Pagell i2004t), building on work by 
iDenicolo et all l|2002fl . iKewlev fc Dopital (j2002D demon- 
strate theoretically how N2 depends on metallicity using 
detailed photoionization modeling. The N2 metallicity 
dependence stems from several effects. As the metallicity 
increases, the ionizing stellar spectrum softens, and the 
overall ionization decreases, as demonstrated by the in- 
creasing ratio of [Oil] /[OIII]. As the ratio of [Oil] /[OIII] 
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Fig. 3. — Two dimensional spectra for objects with spatially resolved Ha emission. Both the spatial scale and positions of Ha and 
[Nil] emission lines are marked. (Left): The spectrum of BX528, which exhibits a tilted Ha emission line with a shear of v max — v m j n ~ 
140 km s — 1 . (Right): The spectrum of BX663, which has a complex Ha emission structure consisting of two components that are offset 
by 390 km s —1 in velocity. The higher-redshift component is spatially extended towards the left of the continuum, while the lower-rcdshift 
component tilts up and to the right of the continuum. This double spatial structure is mirrored in the much fainter [Nil] emission. 



increases, so does the ratio of [NII]/[NIII]. Furthermore, 
at (O/H) metallicities above ~ 0.2 x (0/H) Q , the ratio 
of N/O increases due to the seconda ry nature of Nitro- 
gen in this regime ijHenrv et al.12 000). Limitations of the 
N2 method include the fact that N2 is sensitive to the 
value of the ionization parameter as well as to metallic- 
ity. Also, the N/O ratio can vary as a function age or 
star-formation his t ory (n ot simply abundance). Finally, 
IKewlev fc Dopital l|2002|) show with detailed photoion- 
ization models that while N2 is a monotonically increas- 
ing function of metallicity at less than or equal to solar 
(O/H), it saturates at solar metallicity and even turns 
over towards lower values (at Z > 2Z@). 
Despite these limitations, the N2 indicator offers us 



rough estimates of the H II region abund ances, especially 
using the latest empirical calibration of iPettini fc Page! 
( 2004) . This calibration only includes (O /H) determina- 
tions from the direct T e method or detailed photoioniza- 
tion modeling. The linear relationship between log(0/H) 
and N2 thus derived is: 

12 + log (O/H) = 8.90 + 0.57 x N2 (1) 

For any individual N2 measurement, equation 2] should 
yield the correct (O/H) abundance to within ±0.18 dex 
(68% confidence). Using this relationship, and assum- 
ing its validity for interpreting the integrated spectra of 
high-rcdshift star- forming galaxies, we determine (O/H) 
abundances for the sample of z ~ 2 galaxies with NIR- 
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Table 2. Ha and [Nil] Measurements and the Abundance of Oxygen 



Galaxy 






F Ho c 


*[NII] C 


[Nil]/ Ha 


12 + log(0/H) d 


12 + \og{0/H) c 






( km s ) 








^rCLTlIll oZ, Jr^clgCl J 


(Denicolo) 


Q1623-BX274 


2.4100 


121 ± 10 


9.5 ±0.4 


1.6 ±0.4 


0.17 ±0.04 


8.47 ±0.19 


8.56 ±0.21 


Q1623-MD66 


2.1075 


120 ±5 


19.7 ±0.4 


3.4 ±0.4 


0.17 ±0.02 


8.47 ±0.18 


8.57 ±0.20 


Q1623-BX344 


2.4224 


92 ±9 


17.1 ±0.7 


6.2 ±0.8 


0.36 ±0.05 


8.65 ± 0.18 


8.80 ±0.20 


Q1623-BX453' 


2.1816 


61 ±6 


13.8 ±0.3 


4.1 ±0.3 


0.30 ±0.02 


8.60 ±0.18 


8.74 ±0.20 


Q1623-BX513 


2.2473 


... g 


3.3 ±0.3 


< 0.9 


< 0.30 


< 8.60 


< 8.73 


Q1623-BX528 


2.2682 


142 ± 19 


7.7 ±0.5 


1.5 ±0.5 


0.19 ±0.07 


8.49 ± 0.20 


8.60 ±0.23 


Q1623-BX599 


2.3304 


162 ±9 


18.1 ±0.6 


4.7 ±0.6 


0.26 ± 0.03 


8.57 ±0.18 


8.69 ±0.20 


Q1623-BX663 


2.4333 


132 ± 15 


16.8±0.9 h 


3.5±0.3 h 


0.43 ± 0.05 h 


8.69 ± 0.18 


8.85 ±0.20 



a Vacuum heliocentric redshift of Ha emission. 

b Ha velocity dispersion obtained by fitting a Gaussian profile to the Ha line and deconvolving the effects of instrumental resolution. 
c Line flux and random error in units of 10~ 17 erg s — 1 cm -2 . While the systematic flux uncertainties are ~ 25%, the uncertainty in the 
[Nil] / Ha flux-ratio is determined by the random errors in both line fluxes. 

d Oxygen abundance deduced from the relationship presented in Pcttini & Paecl (2004). For comparison, the most recent estimate of the 
solar abundance is 12 + log(0/H)(n = 8.66 ±0. 05, and that of the Orion nebula is 12 + log(0/H)o rio „ = 8.64 ± 0.06 lAllende Prieto et al.l 
l200llXsp1und et al.ll20oltlEsteban et alJll9981) . 

e Oxygen abundance deduced from the relationship presented in[Dcnicol6 ct al. 1 2 0021) . 

f The Ha and [Nil] line fluxes presented are integrals under Gaussian fits to the lines, whose central wavelengths and FWHM are 
determined by the parameters of the Ha line. For all objects except BX453, the fluxes obtained from the fits to the lines agree with the 
non-parametric integrals under the spectra, well within the uncertainties. However, in the case of BX453, the [Nil] flux obtained from the 
fit is 30% lower than that obtained by integrating non-parametrically under the spectrum, due to the larger apparent FWHM of the [Nil] 
line than the Ha line. 

g The Ha emission from BX513 falls directly on top of a skyline, preventing a measurement of a. 

h The Ha line flux listed for BX663 represents the sum of the two components integrated over the entire extended region of Ha emission. 
The [Nil] line flux and [Nil] / Ha ratio correspond to the more-significantly detected, higher-redshift component, and only include flux from 
the spatial extent common to both transitions. 



SPEC [Nil]/ Ha measurements. The (O/H) abundances 
and associated uncertainties are listed in Table [21 For 
comparison with other recent work, we also list the 
(O/H) abundances im plied by the relation presented in 
iDenicolo et"afl l)2002|) 

12 + log (O/H) = 9.12 + 0.73 x N2 (2) 

which a re systematically highe r than those computed 
from the lPettini fc Paeell l)2004|) relation. In the remain- 
ing discussion, however, w e emphasize the results from 
the iPettini fc Pagell l|2004D calibration, since it is based 
on a more homogeneous sample of (O/H) abundances. 

All of the galaxies for which we detected [Nil] have 
(O/H) abundances consistent with the most recent esti- 
mates of the solar abundance, 12 + log(O/H) = 8.66 ± 
0.05, and tha t of the Orion nebula, 12+log(Q/H)o r .,; rm = 
8.64 ± 0.06 llAllende Prieto et all 120021 TAsrjlund et all 
l2004UEsteban et allll998[) . Furthermore. [Nil]/ Ha sat- 
urates at (O/H) > (0/H) Q and even turns over to- 
wards lower values as the increased cooling at the high- 
est metallicities lowers the H II region T e to the point 
that the collisionally-ex cited [Nil] emission is reduced 
l|Kewlev fc Dopitai r2002'l . Since equation is a linear fit, 
which does not take this turn-over into account, the high- 
est [Nil] / Ha measurements in the current sample - those 
of BX599, BX453, BX344, and BX663 - could very well 
correspond to su per-solar metallicity va lues. To address 
this uncertainty, IPettini fc Paeell l)2004|) propose the use 
of the "03-/V2" indicator, which consists of the ratio of 
two line ratios: ([O III] A5007/ H/?)/([N II] A6583/Ha). 
While [Nil]/ Ha saturates at high metallicity, [OIII]/H/3 
continues to decrease such that 03N2 also continues 
to decrease. The monotonically decreasing behavior of 
03N2 with increasing metallicity has been calibrated 
from 0.3 < Z/Zq < 2.5, and exhibits significantly less 



scatter than the corresponding relation using N2 alone. 
A clear next step is to obtain H-b&nd [OIII]/ H/3 spectra 
for the galaxies in this sample. 

In the preceding discussion, we have interpreted the 
[Nil] / Ha ratios in terms of the H II region oxygen abun- 
dances that they imply. This interpretation rests on 
the assumption that the observed emission lines orig- 
inate in H II regions excited by the photoionizing ra- 
diation from massive O and B stars. The contribution 
from diffuse ionized gas (DIG) to the integrated emission- 
line spectra of local star-forming galaxies gas can en- 
hance the [Nil]/ Ha ratio by different amounts rang- 
ing from less than 0.1 dex to as much as a factor of 
two, r elative to what is measured in individual H II re- 
gions l|Lehnert fc Heckmanlll99l . The DIG component 
could potentially bias our interpretation of [Nil] / Ha ra- 
tios towards higher metallicities. Currently, we have no 
constraints on the contribution of DIG to the spectra 
presented here, which would require spatially-resolved 
Ha narrow-band imaging of the cosmologically-distant 
galaxies in our sample. 

Other possible excitation mechanisms include pho- 
toionization by an AGN power-law continuum, or shock- 
heating, as seen in low-ionization nuclear emission re- 
gions (LINERS) and supernova remnants. In the future, 
combining [Nil]/ Ha with measurements of [OIII]/ H/3 
will clearly discriminate among the different excitation 
mechanisms for these galaxies. For now, there is no ev- 
idence pointing to an AGN as the source of excitation. 
In addition to the NIRSPEC rest-frame optical spectra, 
we have also analyzed the Keck I/LRIS rest-frame UV 
spectra of these galaxies. An AGN manifests itself in 
the rest-frame UV with strong Lya emission, and high- 
ionization emission lines including Si IV, Ci IV, and N V 
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(ISteidel et all20 02). Only one of the seven objects in the 
K s selected sample, BX663, displays strong Lya emission 
in its rest-frame UV spectrum, but no other emission 
lines are detected. In the remainder of the sample, Lya 
is detected in either strong absorption or with an equiv- 
alent width of ~ A, and none of these rest-frame UV 
spectra show evidence for high-ionization emission lines. 

A consistency check on the possible AGN contribu- 
tion is provided by the NIRSPEC spectra. Different 
types of extragalactic emission-line regions can be sep- 
arated according to the relati ve values of the i r emi ssion- 
line strengths. For example , Baldw in et alJ (|1981fl and 
iVeilleux fc Osterbrocld l|1987ft show that AGN and LIN- 
ERS tend to have higher [Nil] / Ha ratios (i.e. on the or- 
der [Nil]/ Ha ~ 1) than star-forming galaxies, and that 
AGN have higher [OIIIJ/H/3 ratios than star- forming 
galaxies for a given [Nil]/ Ha ratio. While we cur- 
rently do not have ii-band [OIIIJ/H/3 spectra for the 
z ~ 2 sample, another way in which AGN and star- 
forming galaxies can be distinguished is according to the 
[01] 6300/ Ha ratio. AGN and shock-heated regions have 
[Ol] 6300/ Ha ~ 0.1 — 1, whereas star-forming galaxies 
have [OI] 6300/ Ha ratios of < 0.1. Since the detection of 
[OI] at less than 0.10 times the strength of Ha pushes the 
limit of what can be achieved with the individual galaxy 
spectra, we construct higher-S/N composite NIRSPEC 
spectra of the entire sample of seven galaxies, and the 
subset of four of these with the highest [Nil] / Ha ratios. 
[OI] 6300 is not detected in either case, yielding a limit of 
[OI] 6300/ Ha < 0.05. The lack of significant [OI] 6300 
emission in the composite NIRSPEC spectra lends fur- 
ther support to the interpretation of the emission line 
strengths in the context of stellar H II regions. 

5. DISCUSSION 

5.1. The Metallicity- Luminosity Relation 

We have presented the rest-frame optical spectra of a 
sample of UV-selected star-forming galaxies that are also 
bright in the rest-frame optical. The [Nil]/ Ha line ra- 
tios of the majority of these galaxies indicate that they 
have already been enriched to at least solar oxygen abun- 
dance at a lookback time of ~ 10 Gyr. It is worthwhile 
to consider their metallicities and luminosities along with 
those of local star-forming and elliptical galaxies. Over a 
factor of 100 in (O/H) abundance and 11 magnitudes in 
blue luminosity, nearby star-forming spiral and irregular 
galaxies exhibit a stro ng trend of hig her metallicity with 
increasing luminosity iJGarn ctt 2002). Elliptical galaxies 
show a similar correlation between metallicity and lumi- 
nosity, although their metallicities are measured from ab- 
sorption indices in l ate-type stars rather th an H II emis- 
sion line strengths (jBrodie fc Huchralll991^ . This corre- 
lation can be explained if less luminous (and less massive) 
galaxies have larger gas fractions, either because they 
have lower star-formation rates per unit mass than more 
massive galaxies, or because they are younger. Another 
way of explaining the correlation is through galactic-scale 
outflows, which could remove a larger fraction of gas and 
metals from less massive galaxies, which have shallower 
gravitational potenti als, thereby r educing the effective 
heavy el ement yield |Garne ttil2002[). 

Bot h IKobulnickv fc Kool pOOqi and IPettini et all 
l)2001|) have shown that z ~ 3 Lyman Break Galaxies 
(LBGs) are ^2 — 4 magnitudes overluminous for their 



metallicities, when compared with the local metallicity- 
luminosity relation. Despite large uncertainties in LBG 
metallicities, the upper bound of the abundances allowed 
by the observed LBG R23 values still falls a factor of three 
in metallicity below the local metallicity-luminosity rela- 
tionship. In Figure0]we plot the l ocal meta l licity- MR re- 
lationship (with data compiled bv lGarnettl l)2002j) ') along 
with the shaded region allowed for LBGs. In this figure, 
we also include the luminosity and metallicity measure- 
ments for the rest-frame optically luminous 2.1 < z < 2.5 
star-forming galaxies presented here. For both the LBGs 
and z ~ 2 galaxies, Mb was calculated using best-fit 
model spectral energy distributions (see section IS~2")l to 
interpolate the observed broadband U n GlZK s colors to 
rest- frame B. 

While t he specific (O/H) value s and errors implied by 
the linear IPettini fc Paeell f 2004) relationship place the 
z ~ 2 galaxies below the metallicity-luminosity relation- 
ship for local galaxies, it must also be borne in mind that 
our estimates of (O/H) represent lower limits because of 
the relative insensitivity of the N2 index when (O/H) 
is higher than solar. Furthermore, since (O/H) abun- 
dances in this plot are derived with the R23 method for 
the z ~ 3 and z = galaxies and with the N2 method 
for the K s < 20.0 z ~ 2 galaxies, it is important to 
acknowledge possible systematic differences between the 
two methods. As mentioned above, these systematic dif- 
ferences are definitely relevant in the super-solar metal- 
licity regime where [Nil] / Ha is no longer sensitive to 
increasing (O/H) abundance. The oxygen abundance 
derived from [Nil]/ Ha may also be ~ 0.1 — 0.2 dex 
lower than the abundance derived from R23 at lower 
metallicities, where [Nil] / Ha is sensitive to metallicity 
l|van Zee et al.lll998j) . To estimate (O/H) from the z ~ 2 
[Nil] Ha ratios in a manner more consistent with the R23 
abundance determinations for local H II regions, we ex- 
amined the emission line data for a set of H I I regions 
located in eight galaxies in the iGarnettJ l)2002j) sample. 
Uniformly measured N2 and R23 values for these H II re- 
gions are contained in Ivan Zee et alJ l|1998|) . The z ~ 2 
galaxies in our sample have [Nil]/ Ha ra tios that span 
from 0.17 to 0.43. H II regions in the Ivan Zee et alJ 
(1998) sample with [Nil]/ Ha in this range are found 
to have 12 + log (O/H) = 8.50 - 9.40, us i ng the upper- 
branch R23 calibration of iZaritskv et al.l l)1994|) . which 
corresponds to 0.7 — 5.5 x (O/H) . While the R23 cal- 
ibration is in itself untested at these very high values 
of metallicity and may well suffer from its own system- 
atic uncertainties, it seems very likely that the z ~ 2 
star-forming galaxies we have observed in this work are 
enriched to levels comparable to those of the most metal- 
rich and massive spiral and elliptical galaxies in the local 
universe. 

Whether solar or super-solar, the metallicities of these 
galaxies are likely to grow further from z ~ 2 to the 
present time, since we observe them while they are 
still in the process of actively forming stars at rates of 
24 — 60Mf7) y r" 1 . In one possible scenario, proposed by 
IMaier et alJ ((2004) to explain the range of (O/H) and 
Mb measurements of star- forming galaxies at z ~ 3, 
z ~ 0.5, and z ~ 0, the rest-frame optically bright z ~ 2 
galaxies in our sample will evolve both in metallicity and 
luminosity by z ~ 0; their oxygen abundances will in- 
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closer to the relationship than z ~ 3 LBGs (blue shaded region). The majority of the z ~ 2 data points represent lower limits in (O/H), 
since the measured N2 values lie in a regime where this line ratio becomes insensitive to increasing metallicity. 



crease by a factor of two, and their luminosities will fade 
by ~ 2 magnitudes. Unfortunately, at this point, we 
cannot precisely determine the amount of future star- 
formation and enrichment in these galaxies (and there- 
fore tracks in metallicity-luminosity space as a function of 
rcdshift) since we have no information about their molec- 
ular gas content or future episodes of gas accretion and 
infall. 

In any case, it is important to obtain deeper [Nil] 
and Ha measurements of K s > 20 star-forming z ~ 2 
galaxies in order to determine the z ~ 2 metallicity- 
luminosity relationship over a larger range of rest-frame 
optical luminosities and relate star-forming z ~ 2 galax- 
ies to those at lower redshift. As a first step, we com- 
p are galaxies in the current sample with those presented 
in lErb et al.l (|2003), in terms of their rest-frame optical 
luminosities (corrected for Ha emission) and [Nil] / Ha 
ratios. The current sample of K s < 20.0 galaxies is char- 
acterized by (M opt ) = -23.29 and ([Nil]/ Ha) = 0.27, 
corresponding to (12 + log (O/H)) = 8.58. Of the 17 
galaxies with H a emission-line flux measurements in 
lErb et aTl l|2003|) . 14 also have measured K s magnitudes, 
from which rest-frame optic al luminosities can be calcu- 
lated. The lErb et alJ |2003) galaxies with K s measure- 
ments are fainter on average than the K s < 20.0 sample, 
with (M op t) = —22.16. It was not possible t o estimate 
[Nil] / Ha on an individual basis for most of the lErb et al.l 
( 2003) galaxies, since [Nil] is not significantly detected 
in the majority of the spectra. To determine the average 



[Nil]/ Ha, we constructed a composite NIRSPEC spec- 
trum from the 14 galaxies with K s measurements. In 
this spectrum, we found [Nil]/ Ha=0.10, corresponding 
to 12 + log (O/H ) = 8.33. Therefore, the galaxies from 
lErb et all {20031 are ~ 1.1 mag fainter in the rest- frame 
optical and 0.25 dex lower in log (O/H) than the current 
sample. This trend marks the preliminary establishment 
of a metallicity-luminosity relationship at z ~ 2. 

5.2. The Masses of Metal-Rich Star-forming z ~ 2 
Galaxies 

The high H II region metallicities of the if s -bright 
star-forming z ~ 2 galaxies suggest mature systems 
that have processed a significant fraction of their bary- 
onic material into stars. Modeling the broad-band op- 
tical and near-infrared (U n G1ZK s ) spectral energy dis- 
tributions of these galaxies with population synthesis 
codes provides an independent test of t his hypothesis. 
Follow ing a si milar procedure to t h at of IShaplev et alJ 
(12001(1 . we use iBruzual fc CharlotJ l)2003|1 models with 
solar-metallicity and a Salpeter initial mass function 
(IMF) to deduce stellar masses a nd ages. Dust e xtinc- 
tion is taken into account with a ICalzetti et al.l (2000) 
starburst attenuation law. 4 In order to determine how 
well the stellar masses and ages can be constrained, 

4 IHeddv fc Steidell 120041) have shown that using the Calzetti law 
to infer unobscured star-formation rates from UV colors and mag- 
nitudes accurately predicts the average X-ray and radio continuum 
fluxes of z ~ 2 star-forming galaxies. 
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Table 3. Star-formation Rates and Stellar Population Parameters 









Uncorrected 


Uncorrected 


Model 














SFR Hq c 


SFK uv d 


SFR C 




Ages 


M t h 
lvl star 


Galaxy 


2Hc a 


Ly. a h 


(M yr-!) 


(MQyr- 1 ) 


(Meyr- 1 ) 


e(b- vy 


(Gyr) 


(lO n M ) 


Q1623-BX274 


2.4100 


4.3 


34 


28 


75 


0.12 


1.3 


1.9 


Q1623-MD66 


2.1075 


6.5 


51 


10 


65 


0.23 


0.9 


0.9 


Q1623-BX344 


2.4224 


7.9 


62 


8 


19 


0.20 


1.6 


1.9 


Q1623-BX453 


2.1816 


4.9 


39 


16 


174 


0.27 


0.4 


0.9 


Q1623-BX528 


2.2682 


3.0 


24 


18 


44 


0.11 


1.7 


1.9 


Q1623-BX599 


2.3304 


7.6 


60 


22 


49 


0.10 


1.3 


1.3 


Q1623-BX663 


2.4333 


7.8 


62 


12 


33 


0.13 


2.0 


2.3 



a Vacuum heliocentric redshift of Ha emission. 

b He* luminosity in units of 10 42 erg s — 1 . ^^^^^^^ 

C SFR calculated from Lh„, using the conversion of Kcnnicutt 1 1998). 

SFR calculated from L150O1 probed by the G apparent magnitude and using the conversion of lKennicutH IT998TV 
C SFR calculated by fitting a Bruzual & Chariot] 120031^ r = 1 Gyr model reddened by dust extinction to the U n G r R.K a colors of the 
galaxies. The best-fit value of E(B — V) allows for the calculation of the dust-corrected value of the star-formation rate, listed here. 
f The best-fit value of E(B — V), assuming a Calzctti ct al. (2000) dust-attenuation law as a function of wavelength. 

g The best-fit value of stellar population age associated with the current episode of star formation, assuming an exponentially declining 
star-formation history with r = 1 Gyr. 
h Formed stellar mass computed by integrating the best-fit exponentially-declining SFR between t = and the best-fit age. 



we investigated a range of star-formation histories of 
the form SFR(t) oc exp(— t/r), with e-folding times of 
r = 0.01, 0.05, 0.10, 0.2, 0.5, 1, 2, and 5 Gyr, as well 
as continuous star- formation models. Before modeling 
the colors, we corrected the K s magnitudes for the ef- 
fects of Ha and [Nil] nebular emission. 5 

Five of the seven galaxies show evidence of having sus- 
tained active star formation over long timescales, are not 
well described by steeply declining star-formation histo- 
ries, and appear to have formed at least ~ 10 11 Mq of 
stars. Only r > 200 Myr models provide acceptable fits 
to their SEDs, yielding ages that are significant fractions 
of the Hubble time at z ~ 2. In fact, while constant- 
star-formation models provide statistically excellent fits 
to the photometry for these five galaxies, the associated 
best-fit ages are older than the ages of the Universe at 
each galaxy's redshift. To avoid this contradiction, we 
only consider models which yield best-fit ages causally al- 
lowed by the finite age of the universe. The r = 200 Myr 
models yield stellar masses of 0.8—1.3 x 10 U M© and ages 
of 0.6 — 0.8 Gyr, whereas the maximally-allowed r models 
for each object yield masses of 1.9 — 2.8 x lO n M0 and 
ages of 2.1 — 2.5 Gyr. Both ground-based J and f/-band, 
and Spitzer rest-frame near- infrared photometry will help 
to discriminate between r = 200 and r > 1 Gyr models. 
The colors of the remaining two galaxies, BX453 and 
MD66, provide very little restriction on star- formation 
history. Both of these galaxies can be fit by steeply- 
declining or constant star- formation models. The best- 
fit masses for BX453 range from 0.4-0.9 x 10 n Af o with 
associated ages ranging from 0.1 — 0.5 Gyr, whereas the 
masses for MD66 range from 0.5 — 1.2 x IO^Mq and ages 
from 0.3 — 1.7 Gyr. While less massive than the majority 
of the sample, BX453 and MD66 have already formed a 
significant fraction of the stellar ma ss of an L* galaxy in 
the local universe llCole et alJl200l . 

For all seven galaxies, the uncertainties in star- 
formation history lead to systematic uncertainties in the 

5 For a full de scription of the modeling procedure, see 
IShaplev et"aH <200ll) . 



stellar mass of a factor of ~ 2. Furthermore, if there is 
an underlying maximally-old population, hidden at rest- 
frame UV-to-optical wavelengths by the current episode 
of star formation, the stellar masses derived here could 
underestima te the total stellar m ass by up to a fac- 
tor of five ijPapovich et alJ l2001|) . Both the inferred 
dust extinction, E(B — V), and dust-corrected star- 
formation rate depend systematically on the assumed 
star-formation histories; specifically, the best-fit mod- 
els are characterized by more dust-extinction and higher 
bolometric star- format ion rates when we assume larger 
t values. Assuming the maximally-allowed r for each 
galaxy yields E(B - V) = 0.12 - 0.27 for the sample, and 
SFR = 30 - 180Af Q yr" 1 . In Table we list the best- 
fit stellar-population parameters, assuming t = 1 Gyr, 
which provides causally allowed ages for all galaxies in 
the sample. FigureElshows plots of the best-fit r = 1 Gyr 
models. While these objects were selected to be at z ~ 2 
because they display the colors of star- forming galaxies, 
it is clear that they are not only metal-rich, but also have 
large stellar masses. 

The z ~ 2 star -forming galaxies in the sample of 
iSteidel et all l|200l were selected to have intrinsic rest- 
frame UV colors very similar to those of star-forming 
LBGs at z ~ 3. In relating the two samples of galax- 
ies, it is important to consider if plausible progenitors of 
the massive, metal-rich K s < 20.0 "BX" galaxies can be 
found among the sample of z ~ 3 LBGs with near-IR 
photometry, whose s tellar populations were modeled by 
IShaplev et al.1 II2001D. T o addr ess this question, we re- 
model the IShaplev et pOOl sample, using the most 
current relea se of the Bruzual fc Chari ot population syn- 
thesis codes ijBruzual fc Charlotil2003j) . For a consistent 
comparison with the best-fit values presented in Table El 
we assume exponentially declining star-formation histo- 
ries, with t = 1 Gyr. Evolving all the z ~ 3 best-fit 
models forward in time, we compute the star-formation 
rates and stellar masses formed by z = 2. Assuming that 
E(B — V) does not change between z ~ 3 and z — 2, we 
also compute the predicted 1Z and K s magnitudes, and 
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Fig. 5. — Best-fit Bruzual & Chariot (2003) models with r = 1 Gyr. These plots correspond to the stellar population parameters listed 
in Table|3 U n GRK B data points are plotted, as are photometric error bars, which are typically smaller than the size of the data points. 
The best-fit models are all characterized by moderate amounts of dust extinction and significant Balmer Breaks, indicating the extended 
time-scale over which star-formation has proceeded and 1O 11 M0 of stellar mass has built up. 



U n — G, G — 1Z colors. The predicted UV colors can 
be used to determine if a galaxy would be r ecover ed by 
the z ~ 2 UV criteria of lAdelberger et all l|2j)0_4f), and 
the predicted near-IR magnitudes indicate whether or 
not a galaxy would be included in the rest-frame opti- 
cally bright sample presented here , with K K < 20 .0. At 
least 15% of the z ~ 3 LBGs fit bv lShaplev et~al1 (Effllh 
have stellar population models time-evolved to z — 2 
with stellar masses, star-formation rates, E(B — V) and 
age within the range of properties of the K s < 20.0 BX 
galaxies presented here. The time-evol ved U n GTZ colors 
of thes e LBGs satisfy the UV criteria of lAdelberger et aD 
()2004[K and they also have K s < 20.0. Therefore, at 
least 1 5% of the z ~ 3 LBGs modeled bv lShaplev et al.l 
l)2001() constitute plausible progenitors for the fsT s -bright 
BX objects. As shown below, the if s -bright BX's have 
a comoving number density that's only 14% as high as 



that of LBGs in the iShaulev et alJ l)2001[) sample, so we 
conclude that most of them could have been LBGs in the 
past. 

It is now possible to estimate the stellar mass density 
associated with K s < 20.0 UV-selected z ~ 2 objects. 
First we estimate the average mass of the 7 objects with 
[Nil]/ Ha measurements. Taking into account the sys- 
tematic uncertainties due to the range of allowed star- 
formation histories, we find (M st ar) — 1.4±0.5x 1O 11 M . 
This average mass should be representative of the entire 
sample of K s < 20.0 z ~ 2 galaxies since the distribu- 
tions of G — 1Z and 1Z — K s colors of the seven galaxies 
with mass estimates are similar to those of the entire 
sample of rest-frame optically luminous objects. 

Next we estimate the comoving space density of K s < 
20.0 z - 2 UV-selected galaxies. In the Q1623 WIRC 
pointing, which covers 92 arcmin 2 , there are 87 "BX" 
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Table 4. Fraction of 2.0 < z < 2.7 BX objects vs. TL magnitude 

TZ Magnitude Range Fraction of BX objects with 2.0 < z < 2.7 a Q1623 BX objects with K s < 20.0 and no z 

Tl < 21. 7 b 0/48 39 

21.7 <1l< 23.0 6/51 13 

23.0 < n < 24.0 145/301 7 

24.0 < Tl < 25.5 386/556 8 

a These statistics are based on the entire sample of 956 spectroscopically identified BX objects. 

b We use TZ = 21.7 as a cut-off since all spectroscopically-confirmed BX objects at brighter magnitudes are stars and QSOs. 



objects with K s < 20.0, including 13 galaxies with con- 
firmed redshifts between 2.0 < z < 2.7, one broad-lined 
QSO at z = 2.53, four stars, and two low-redshift galax- 
ies at z ~ 0.2. The remaining 67 photometric z ~ 2 
candidates have either not been attempted spectroscopi- 
cally or have been attempted unsuccessfully, with unde- 
termined redshifts. To estimate the number of K s < 20.0 
galaxies without spectroscopic redshifts that are likely to 
fall in the redshift range of interest, we use the redshift 
statistics as a function of 1Z magnitude from our total 
samp le of 956 BX obje cts with spectroscopic informa- 
tion (jSteidel et alJ l20041. Based on the BX high-redshift 
fractions and the distribution of 1Z magnitudes for the 67 
K s < 20.0 sources without redshifts presented in Tabled] 
(i.e. multiplying the second column of Table 0] by the 
third column), we estimate that an additional 11 K s < 
20.0 galaxies without measured redshifts are likely to be 
at 2.0 < z < 2.7. The total number of K s < 20.0 BXs at 
2.0 < z < 2.7 is therefore N = 13 + 11 = 24. Weighted . 
by the BX redshift selec tion function (jSteidel et al.l200"4l 
lAdelberger et alJl2004|) . the comoving volume between 
2.0 < z < 2.7 is 1496 Mpc 3 arcmkr 2 . Over the en- 
tire Q1623 WIRC pointing, the space-density of 2.0 < 
z < 2.7 UV-selected objects with K s < 20.0 is therefore 
n = 1.7 x 10~ 4 Mpc~ 3 . The space density of massive 
galaxies is of course prone to large sample variance be- 
cause such objects are expected to be strongly clustered. 
While not yet properly quantified, initial statistics from 
our two other fields with i^T s -band data suggest variations 
of a factor of ~ 2 from field to field. 

Multiplying this space density by the average mass of 
the rest-frame optically luminous objects we obtain p = 
n x (M star ) = 2.5 ± 0.9 x 1O 7 M Mpc" 3 , whi ch is 4% of 
the z = stellar mass density determined bv lCole et alJ 
(|200H) . The K s < 20.0 objects represent only ~ 10% of 
the sample of BX objects at z ~ 2, and since we expect 
to find massive systems down to significantly fainter K s 
magnitudes, this figure represents an extreme lower limit 
on the stellar mass density associated with UV-selected 
star- forming galaxies at z ~ 2. In the future, when we 
have the distribution of stellar masses over a much wider 
range in luminosity in hand, we will determine a more 
robust constraint on the stellar mass densit y at z ~ 2, 
extending the work of iDickinson et aD l|200 3) to a larger 
cosmic volume. 

5.3. Comparison with Other Surveys 

There are other current surveys targeting massive 
galaxies in the same redshift range. It is instructive to 
compare the properties of the K s < 20.0 UV-selected 
objects presented here with those of the galaxies in the 



K20 survey ijCimatti etjtljjgj, ) Q 2|) and the G emini Deep 
Deep Survey CODDS: lAhraham et al.ll200l . The K20 

project is a spectroscopic survey of ~ 500 K s < 20.0 ob- 
jects, 9 of which have been spectroscopical ly confirmed 
in the redshift range 1.7 < z < 2.25. iDaddi et all 
(2004) demonstrate that this subsample of galaxies is 
characterized by dust-corrected star-formation rates of 
~ 100 — 500M Q yr -1 , constant star-formation ages of 
0.25 - 1.7 Gyr, and stellar masse s of 0.3 - 5.5 x 1O U M . 
Six of the galaxies in lDaddi et al.l l)2004|) are at z > 2, and 
four of them have op tical colors that s a tisfy the "BX" 
criteria presented in lAdelberger et all l)2004j) . These 
four galaxies have E(B — V) = 0.3 and constant star- 
formation ages of 0.7 — 1.7 Gyr, very similar to the 
properties of BX453 and MD66 when constant star- 
for mation histo r ies ar e assumed. Of the three galaxies 
in IDaddi et alJ l)2004|) that are at z < 2, two would- 
be rec overed by the "BM" criteria of lAdelberger et "all 
(2004), which target star-forming 1.5 < z < 2.0 galaxies. 
Since the K20 s urvey uses a differen t optical filter set 
from that of the iSteidel et alJ (I2004D survey, we calcu- 
lated the U„GTZ colors of the IDaddi et alJ l|2004fl galax- 
ies by passing the best-fit redshifted, dust-reddened, con- 
stant star-formation models (from their Table 1) through 
the U n GlZ filter set. There appears to be significant 
overlap between the K s < 20.0 UV-selected objects 
and the high-redshift tail of the K20 survey. Further- 
more, the space density of 1.7 < z < 2.25 K20 objects 
is 1.6l g x 10 -4 Mpc -3 , which is very comparable to 
the space density of K s < 20.0 UV-selected objects at 
2.0 < z < 2.7 in the Q1623 field, n = 1.7 X 10~ 4 Mpc" 3 . 
Finally, the distribution of stellar masses of the K20 and 
rest-frame optically luminous UV-selected objects are 
very similar (with the exception of one K20 object that 
has 5.5 x lO n M0 and would not be selected as a BX). 
The additional UV-selection criteria appear to recover 
the majority of K20 objects, because most > lO 11 ^/© 
galaxies harbor active star formation at z ~ 2. 

While the majority of the 300 galaxies in the GDDS 
survey are at z < 1.5, this infrared-selected K < 20.6 sur- 
vey cont ains a small high-redshi ft tail of 1.6 < z < 2.0 
galaxies l)Glazebrook et al J 120041) . Using GDDS galax- 
ies to tra ce the evolution of stell ar mass-density in the 
universe. IGlazebrook eTaTl 1I2OOI find a 1.6 < z < 2.0 
mass density for M > Mu m = 6.3 x 1O 1O M objects of 
3.3 ± 1.2 x 10 7 M© Mpc" 3 . Below M Hm , the K < 20.6 
GDDS becomes incomplete for objects with the highest 
mass-to-light ratios. The GDDS masses and mass den- 
sities were ca lculated assuming an IMF that turns ove r 
below O.5M l|Baldrv fc Glazebrookl2 003: Kro upal200l"|) : 
it is therefore necessary to multiply the GDDS numbers 
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by a factor of 1.82 for a meaningful comparison with 
the results presented here, derived using a Salpeter IMF. 
The range of masses of the spectroscopically confirmed 
GDDS 1.6 < z < 2.0 galaxies is very similar to that of the 
K s < 20.0 UV-selected objects, while the mass density, 
computed down to a limiting magnitude of K < 20.6 
is ~ 3 times higher. A more fair comparison in mass 
densities will result when we have determined the stellar 
masses of fainter UV-selected objects. Nonetheless, it is 
still currently possible to determine how complementary 
the UV-selected and GDDS surveys are in the z > 1.6 
redshift range, since one of the GDDS survey areas, the 
SSA22 field, has also been image d in the U n GTZ filte r set 
as part of the z ~ 3 LBG survey ijSteidel et al.ll2003|) . In 
the SSA22a field, there are three spectroscopically con- 
firmed GDDS galaxies with 1.6 < z < 2. 0, and four with 
uncer tain redshifts of 2.0 < z < 2.2 ijAbraham et al.l 
|2004|) . All three of the 1.6 < z < 2.0 galaxies are re- 
covered with the BX/BM J7 n G72.-selection criteria, and 
three of the four putative z > 2 objects would be selected 
as BXs. At 1.4 < z < 1.6, there are 12 GDDS objects in 
the SSA22a field, half of which would be recovered by the 
BX/BM UV-selection criteria. The 50% overlap between 
GDDS galaxies at 1.4 < z < 1.6 and BX/BM objects is 
not especially meaningful since the BM color criteria are 
intended to select galaxies at 1.5 < z < 2.0, and the 
BX color criteria target galaxies at 2.0 < z < 2.5-i.e. 
1.4 < z < 1.6 lies at the edge of the BX/BM selection 
function, where the selection efficiency drops off. How- 
ever, the large overlap of the GDDS and BX/BM surveys 
at z > 1.6 is again consistent with the idea that even K s - 
band selected surveys that are supposed to be "mass- 
selected" include a significant fraction of star-forming 
galaxies at z ~ 2. 

5.4. Stellar vs. Dynamical Masses 

The preceding discussion relies on population synthe- 
sis models of the galaxy colors to infer stellar masses 
and therefore the comoving stellar mass density in large 
galaxies at z ~ 2. Rest- frame optical spectroscopy al- 
lows for a completely independent probe of galaxy mass: 
the Ha velocity dispersion. A physical scale-length and 
assumption of the geometric configuration are necessary 
to convert the velocity dispersion into a dynamical mass. 
In previous determination s of high-redshift star-forming 
galaxy dyn amical masses llPettini et al.ll200ll lErb et aLl 
120031 l2004j) . a spherically-symmetric distribution of mat- 
ter was assumed, and the physical scale-length adopted 
was the typical z ~ 3 or z ~ 2 galaxy continuum half- 
light radius, as determined from Hubble Space Telescope 
images, r 1 / 2 =0'.'2-0'.'3. From the excellent seeing condi- 
tions during our NIRSPEC exposures, we also determine 
Ha half-light radii of 0'.'2-0'.'3. Using these radii, and 
assuming a spherical geometry, we determine dynamical 
masses with the formula: Md yn — 5cr 2 r 1 / 2 /G. Dynami- 
cal mass estimates are listed in Table as well as stellar 
masses inferred from the smallest t allowed by the galaxy 
photometry, and the largest r allowed by the age of the 
universe constraint. As discussed in section I5~21 the stel- 
lar mass of the best-fit model is an increasing function of 
the assumed time-constant, r. The actual stellar mass 
for each object lies somewhere between the limits listed in 
Tabled (neglecting the mass of an underlying maximally- 
old burst). 



For every galaxy, if a Salpeter IMF extending down 
to O.1M is assumed, the allowed range of stellar mass 
exceeds the calculated dynamical mass within ri/2. For 
BX663, MD66, BX528, and BX599, this discrepancy is 
roughly a factor of ~ 2 - 4, whereas for BX344, BX274, 
and BX453, the stellar mass estimate exceeds the dy- 
namical mass by an order of magnitude. Perhaps this 
discrepancy is not surprising, given that the galaxy stel- 
lar mass was determined from luminosities and colors ex- 
tending over a physical region in some cases several times 
as large in radius as the region enclosed by the half-light 
radius. A larger radius may be more suitable for infer- 
ring the dynamical mass. Indeed, with the assumption 
of a constant M/L ratio as a function of radius, the mass 
enclosed within 7*1/2 should be multiplied by 2 in order 
to compare with the total stellar mass. 

Furthermore, there are observational effects that may 
apply if the geometry is not the simple spherically- 
symmetric one assumed for the calculation. If the emis- 
sion is extended along a specific axis, the possible mis- 
alignment of the slit with this axis must be taken into 
account. If the configuration of H II regions is more 
disk-like, then inclination effects may also cause us to 
underestimate the t he circula r veloc ity, and therefore, 
velocity dispersion. lErb et all l)2003|) estimate that the 
combination of these two effects may result in an un- 
derestimate of a factor of ~ 2.5 in circular velocity. 
Of course, the numerical coefficient multiplying a 2 r de- 
pends on the assumed geometry as well, and is unity 
for a disk geometry (as opposed to 5, assumed for the 
spherical configuration). Independent of slit misalign- 
ment and disk inclination effects, the emission within 
the Ha half-light radius may not sample the full range 
of rotational velocities of the underlying disk. Using local 
starburst galaxies wi t h wel l-determined rotation curves, 
ILehnert fc Heckmanl l|1996f) have shown that the widths 
of optical emission-lines originating in the bright central 
regions of local starburst galaxies only sample the solid- 
body part of the rotation curve, and therefore do not 
reflect the full circular velocity attained o n the flat part 
of the rota tion curve at larg er radii. Both lPettini et aT] 
l)2001|) and lErb et alJ 0003) have pointed out this lim- 
itation in interpreting the rest-frame optical lincwidths 
of z ~ 3 and z ~ 2 galaxies in terms of masses. Finally, 
stellar masses were determined assuming a Salpeter IMF 
extending down to O.IMq, though we have no observa- 
tional constraints on the low-mass regime of the IMF at 
high redshift. An IMF which turns ov er below 0.5M Q , 
as assumed bv lGlazebrook et all l|2004|) . would lower the 
stellar mass estimates by a factor of ~ 2. We assumed 
the standard Salpeter form in section 15.21 for the pur- 
poses of comparing with earlier studies of the stellar mass 
content of galaxies at low and high redshift l|Cole et alJ 
120011 IShaplev et al.ll200lt IPapovich et alJl2001|) . though 
a more realistic form for the low-mass end of the IMF 
would significantly reduce the systematic discrepancy be- 
tween stellar and dynamical masses (as shown by the 
values in parentheses in Table EJ). 

Most of the above systematic effects would have caused 
us either to underestimate the dynamical mass, or over- 
estimate the stellar mass. There is enough uncertainty 
in our estimates of the low-mass end of the IMF, the 
appropriate radius at which to evaluate Md yn , and the 
geometrical configuration of H II regions, that the cur- 
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Table 5. Comparison of Dynamical and Stellar Masses 







Mj b 
lvl dyn 


' mm 


M**„Ut ■ 1 d 

1 * 1 star \ i rmn } 


T C 

' max 


A/f i _ _, /" T\™, ™ ^ 

ilJ scar^' max / 


Galaxy 


(km s" 1 ) 


(lO n M ) 


(Gyr) 


(lO n M ) 


(Gyr) 


(10 n M Q ) 


Q1623-BX274 


121 


0.28 


0.20 


1.27 (0.70) 


5.00 


2.76 (1.52) 


Q1623-MD66 


120 


0.28 


0.05 


0.49 (0.27) 


oo 


1.19 (0.65) 


Q1623-BX344 


92 


0.16 


0.20 


1.06 (0.58) 


2.00 


2.18 (1.20) 


Q1623-BX453 


61 


0.07 


0.01 


0.49 (0.27) 


oo 


0.92 (0.51) 


Q1623-BX528 


142 


0.58 


0.20 


1.09 (0.60) 


2.00 


2.47 (1.36) 


Q1623-BX599 


162 


0.50 


0.20 


0.83 (0.46) 


5.00 


1.69 (0.93) 


Q1623-BX663 


132 


0.45 


0.20 


1.06 (0.58) 


1.00 


2.35 (1.29) 



a Ha velocity dispersion. 

b Mass calculated from the Ha velocity dispersion. 

c Minimum time-constant of models of the form SFR(t) oc exp(— t/r) that provide statistically acceptable fits to the galaxy colors. 

d Mass of best-fit model, assuming r = T m j n . Values not i n par entheses assume a Salpeter IMF extending down to O.IMq. Values in 
parentheses assume the more realistic lBaldrv Glazcbrook (200ij) IMF, which has a break at IMq. 

"Maximum time-constant of models for which the best-fit age is younger than the age of the universe. BX453 and MD66 can be fit by 
r = oo, i.e. constant star-formation models, while the remaining galaxies' best-fit ages exceed that of the universe when r = oo. 

f Mass of best-fit model, assuming r = T max . Values with and without parentheses have the same meaning as defined in note (d). 



rent discrepancy between stellar and dynamical mass es- 
ti mates does not repre sent a crisis. However, the method 
of IPettini et"all ({2001^ and lErb et all <|2003j) . using sim- 
ple assumptions to compute dynamical masses, clearly 
does not account for the amount of stellar mass derived 
from modeling the colors of these galaxies. We must un- 
derstand the nature of these discrepancies before we can 
trace the evolution of mass in galaxies as a function of 
redshift. 

6. SUMMARY 

We have presented the rest-frame optical spectra of a 
sample of seven UV-selected z ~ 2 galaxies that also have 
K s < 20.0. We reach the following principal conclusions: 

1. The majority of these galaxies exhibit high 
[Nil]/ Ha ratios indicative of solar, and possibly super- 
solar, metallicities, at a l ookback time of ~ 10.5 Gyr. 

2. The calibration of IPettini k. Paeell (|2004|) . which 
converts [Nil]/ Ha ratios into (O/H) abundances, indi- 
cates that the galaxies in the sample are over-luminous 
for their metallicities relative to the local metallicity- 
luminosity relationship. However, we may have underes- 
timated the (O/H) abundances due to the insensitivity of 
the A^2 indicator at super-solar metallicity. In any case, 
given the active star-formation rates inferred from Ha lu- 
minosities and population-synthesis modeling, these ob- 
jects will likely become more metal-rich at later times, 
and are probably the progenitors of massive spiral and 
elliptical galaxies in the l ocal universe. T he sample of 
z ~ 2 galaxies presented in lErb et al.l l)2003() has an aver- 
age rest-frame optical luminosity ~ 1.1 mag fainter and 
average [Nil]/ Ha ratio three times smaller, relative to 
the corresponding values for the current sample, consis- 
tent with the existence of a metallicity-luminosity rela- 
tionship at z ~ 2. 

3. Population-synthesis models of their broad- 
band rest-UV/optical spectral energy distributions sug- 



gest that these galaxies have been forming stars over 
timescales comparable to a Hubble time, contain > 
10 u Mq of stars, and are still actively forming stars. 

4. Comparison with surveys such as the K20 and 
GDDS projects indicates a significant overlap between 
near-infrared- and UV-selected objects. A large fraction 
of the most massive galaxies at z ~ 2 are actively forming 
stars. 

5. Dynamical mass estimates derived from Ha veloc- 
ity dispersions and galaxy half-light radii are systemati- 
cally lower than stellar mass estimates deduced from the 
broad-band SEDs of the galaxies presented here. While 
several plausible effects can be invoked to explain the dif- 
ference in mass estimates, this discrepancy must be un- 
derstood if we are to trace the evolution of galaxy masses 
in the universe. 
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